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Abstract. The reconstruction procedure, which has proven quite useful to obtain viable 
models of the universe evolution, is here employed in order to construct inflation models. 
It has the advantages that it ensures full consistency with astronomical observations and 
that it allows to evaluate the stability of the resulting cosmological model quite easily. The 
reconstruction for two different types of Lagrangian, included in the frame of G-inflation, is 
carried out in detail and explicit models for each Lagrangian are constructed. As a bonus for 
having used this reconstruction formalism, the final models are easily adjusted to satisfy the 
observational constraints—imposed by the most recent data releases of the Planck mission— 
on the spectral index, the tensor to scalar ratio, and the running of the spectral index. 
Further, it turns also to be not difficult to make the models stable. Thus, the method here 
developed provides a general and very efficient tool, a quite natural procedure to construct 
models consistent with very precise observations. It can also be applied to other models, 
besides the ones here considered. 
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1 Introduction 

Inflation in the early universe [1-3] has become a part of the standard cosmological model. 
As is known, its role is fundamental for several important reasons, starting with its role in the 
amplification of the incipient quantum density or curvature fluctuations. One problem, as of 
now, is that there are too many models of inflation and one should try to establish which is 
the one that better describes the very first stages of the evolution of our universe. Actually, 
the true nature of inflation is becoming more and more clear thanks to the latest data 
releases corresponding to the very precise observations of the cosmic microwave background 
(CMB) carried out recently by the Planck satellite [4, 5], the BICEP2 experiment [6] and the 
Wilkinson Microwave anisotropy probe (WMAP) [7]. The most common mechanism driving 
inflation involves a scalar field and relies on the specific form of the potential of this scalar 
field, called the inflaton [2, 3, 8]. 

Since the appearance of this potential-driven inflation, a good number of novel models 
have been devised by modifying the kinetic term of the theory. These include k-inflation 
[9], ghost inflation [10], Dirac-Born-Infeld inflation [11], the simplest higher-derivative theory 
without ghosts, which is modified Gauss-Bonnet gravity [12], G-inflation [15], and generalized 
G-inflation [16]. Some of these models are already on the verge of being excluded by the 
observational constraints obtained from the Planck and WMAP data releases, but many 
of them are still compatible with these astronomical observations and thus one cannot yet 
decide which is the preferred inflationary model. 

As an alternative approach in order to build inflation models, the Hamilton-Jacobi 
formalism of inflation has been put forward. It also deals with a canonical scalar field 
model with a potential [17] but in this formalism one provides the Hubble parameter as a 
function of the scalar field by hand, and one can then determine the potential term from 
it. This formalism has been extended to models with a non-canonical kinetic term [18, 19]. 
In parallel with this development, but now mainly in the context of the present acceleration 
of the universe, a formalism which is quite similar to this one, termed the reconstruction 
method, was pioneered in [20] and subsequently developed with great success (for a review of 
cosmological reconstruction, see [21, 22]). In the reconstruction formalism, using re-scaling of 
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the scalar field, one can eventually provide the Hubble parameter as a function of time instead 
of the scalar field, so that it becomes easier to compare the model with the observations than 
in the Hamilton-Jacobi formalism. So far, however, the reconstruction formalism has been 
used in this context in a limited number of models only [23-25]. 

The purpose of the present work is to explicitly develop and carry to the end the 
reconstruction formalism for two different types of Lagrangian, included in the frame of 
G-inflation, for which the reconstruction formalism has not been thoroughly devised. The 
specific form in one of them consists of a generic function of the kinetic term and a potential 
term of the scalar field. The other one exhibits the same form as G-inflation [15]. Once 
we finish with the construction of concrete models for each Lagrangian, we will provide a 
proper Hubble parameter together with a quantity related to the sound speed of the scalar 
field and demonstrate that they satisfy the stringent observational constraints imposed on 
the spectral index, on the tensor to scalar ratio, and on the running of the spectral index, 
by the most recent data releases of the Planck mission [4]. Further, as we will see in what 
follows, it will also be not difficult to make the models stable. We will see how the method 
here developed provides a general, efficient and quite natural procedure to construct models 
which are consistent with the most precise observations. We will do that by working explicitly 
with some examples but the method can also be applied to other models besides the ones 
here considered. 

The paper is organized as follows. In Sec. 2, we give a brief review of Galileon cosmology. 
In Sec. 3, we describe the reconstruction formalism for the two types of Lagrangian considered 
and construct specific models for each of the Lagrangians. Sec. 4 is devoted to conclusions. 
Finally, in an Appendix, we describe the reconstruction formalism for the case of exact de 
Sitter background solutions in G-inflation. 

2 Short review of Galileon cosmology 

In this section we provied a summary of the Galileon theory [13, 14] and its cosmological 
consequences following the notations in [16] . We actually consider only a part of the Galileon 
action, namely 

s = j d 4 xV=g R + Ktf, X ) - G(X)g^X t y v c^j , ( 2 . 1 ) 

where Mp\ = 1.2 x 10 19 GeV is the Planck mass, R the scalar curvature, K and G are generic 
functions of a single scalar field, 0, and X = — ;,</>. We restrict the function G 

to depend on X , only, and fix other Galileon terms, as G 4 = Mpj/2 and G5 = 0, just for 
simplicity. 

For describing the background evolution of the isotropic homogeneous universe, we 
assume that 4> = </>(£) and work with the flat FLRW metric ds 2 = — dt 2 + a(t) 2 dx 2 , where 

a{t) is the scale factor. The Einstein equations and the field equation derived from the action 

(2.1), in this background, are 

3 M^H 2 = 2 XK x -K + 6 X<j)HG x , ( 2 . 2 ) 

Mpj(3iL 2 + 2 H) = -K + 2Xj>G x , (2.3) 

(<j>K x + 6 HXGx)' + 3 H(j>K x + GHXGx) = , (2.4) 
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where H = a/a is the Hubble parameter, the dot meaning time derivative. The indexes on 
the functions K and G denote corresponding partial derivatives. These equations determine 
the background evolution of the universe. 

Next, we will obtain the quadratic action for the scalar and tensor modes of the cosmo¬ 
logical perturbations of the full action (2.1). We fix the unitary gauge <= cj>(t) and perturb 
the metric using the ADM formalism, in the form 

ds 2 = -A f 2 dt 2 + 7 ij -(da; i + N i dt){dx j + N j dt ) , (2.5) 


where 

Af = 1 + a, Ni = di/3, jij = a 2 (f)e 2C fs^ + Kj + ^h ik h k j) . (2.6) 

In the above, a and /5 are the auxiliary fields for the scalar perturbations, £ for the curvature 
perturbations and hij for the tensor perturbations, which satisfy dihij = 0 and ha = 0. After 
removing the auxiliary fields by solving the constraint equations, we obtain the following 
quadratic actions for the tensor and for the curvature perturbations, respectively, 


where 


4 2) = f dtdx 3 a 3 


• 2 _ (Vhjj) 
ij a 2 


si 2) = / dtdx 6 a 


3„3 


QsC 2 - 


(VC ) 2 


?S - M Pl “ M P\ > 

Qs = ^piq 2 + 3Afpj , 

0 = -j)XG x + Mp { H , 

S = XK x + 2 X 2 K X x + 12 H<j)XG x + 6H^X 2 G X x ~ 3 M^H 2 . 


(2.7) 


( 2 . 8 ) 

(2.9) 


The square of the sound speed of the tensor perturbations is the same as the light speed, 
and for the case of the curvature perturbations, we have 


c 


2 

s 


Fs_ 
Qs ' 


( 2 . 10 ) 


We need Xs > 0 and Qs > 0 in order to stabilize the system against the ghost and the gradient 
instabilities. Furthermore, to prohibit superluminal propagation, c 2 > 1, what is compulsory 
for avoiding the closed time-path problem, we require Qs — J~s > 0. To summarize, for both 
stability and non-superluminal propagation, we need to impose 


Xs > 0 and Qs ~ Xs > 0 


( 2 . 11 ) 


To obtain the tensor and the scalar power spectrums, we assume that the following 
parameters are nearly constant: 

e =-^ ~ const, f s = —4 — const, q s = ~ const . (2.12) 

H 2 ’ J HX S HQ S 
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Here, we further assume that all parameters (2.12) are much smaller than 1, namely e, f s , g s <C 
1 1 . After deriving the mode functions from the actions (2.7), the tensor and scalar spectra 
are obtained, up to first order in the parameters (2.12), as 


V h = 




2 H 2 

M 2 ^ 


k=aH 


H 2 


8ir 2 c s F s 


c s k=aH 


(2.13) 

(2.14) 


where A: is a co-moving wave number of the perturbations. It is apparent that the tensor 
spectrum depends on time only through the Hubble parameter. On the other hand, the 
scalar spectrum depends on time not only through the Hubble parameter but also through 
the sound speed of the scalar perturbations and the quantity Fs- The spectral index, up to 
first order in the parameters ( 2 . 12 ), is 

n s - 1 = -2e - ^f s + g s , (2.15) 

and the tensor to scalar ratio reads 

r = 16 irr Cs • ( 2 - 16 ) 

1V1 PI 

The running of the spectral index is provided by 

^2e+ ^f s - g s ^j . (2.17) 

The observational values for these quantities, at the pivot scale k * = O.OSMpc^ 1 , reported 
by the Planck mission [4] are 

ln(10 10 P c ) = 3.062 ± 0.029 , 
ns = 0.9690 ± 0.0063 , 
r < 0.176 , 

a s = —0.0076^QQQgQ, (Planck TT+lowP+lensing) . (2-18) 

Let us now imagine that we have obtained an explicit function of time for H and a quantity 
related to c s from the astronomical observations data, for example, from (2.13) and (2.14). 
An action which will produce both this function H and the quantity related to c s will be 
obviously consistent with the observation. The reconstruction method enables us to construct 
such an action. Using this method one is able in fact to provide an explicit function of time 
for H and the quantity related to c s , which are consistent with observations, and the stability 
condition ( 2 . 11 ) by hand , as inputs of the model we consider, and by using them one can 
derive an action which will produce those inputs as the result of the evolution of the system. 
In the next section, we will formulate the reconstruction procedure for two different types of 
Lagrangian. 

1 In [16], the power spectra of the perturbations and the spectral index were derived without the assumption 

that t,f s ,g a -C 1. Our expressions for the power spectra and the spectral index do match those obtained in 

[16], after imposing these conditions to the corresponding expressions. 
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3 Reconstruction procedure 


In this section, we formulate the reconstruction of the action for two different types of La- 
grangian with an explicit function of H and a quantity related to c s . 

To start, we recall the e-folding number definition, 


N = r Hdt = In— , 

Jt a 

dt dN ’ 


(3.1) 


where t e and a e are the cosmic time and the scale factor when N = 0, respectively. We use 
N as affine parameter, instead of t, because we will construct models of inflation, and in such 
models IV is a more useful parameter to directly constrain model parameters by observations, 
than is t. We assume that inflation ended at N = 0. At the same time, we redefine the scalar 
field as 


0 = 



H(<p) 


dip , 


X = ~\= 


g(^) 


g^dpipdvip , 


(3.2) 


and we assume that we can take ip = IV as one of the background solutions of the system 
just by choosing a field redefinition (3.2) properly. We note that the function q(ip = N) must 
be positive, since we consider a real scalar field theory. In this case, at the background level, 
we get following relations 


/ 


1 - 


VMN) 

H(N ) 


dN , 


X = -j 2 = q(N) , = ^%(N) . 


(3.3) 


The background equations (2.2) using the above relations (3.3) become 


3 M^H(N) 2 = 2qK q (N, q ) - K(N, q) + 6 H^G q (q) , 

M| 1 (3R(AT) 2 - 2 HH\N)) = -K(N,q ) - Hq'^G 3q (q) , 

(y/2qK q (N, q) + 6 HqG g (q)y - 3(^K q (N, q) + 6HqG q (q )) = A ^ ,g) , (3.4) 

where the prime denotes a total derivative with respect to N, and the functions with the 
indices N or q mean partial derivatives with respect to those indices. Next, we will formulate 
the reconstruction explicitly for two different types of Lagrangian. 


3.1 P(X)-V(<j>) 

Consider the following Lagrangian 


K(<i>,X) = P(X) - V{y{ct>)) , G(X) = 0 . 


( 3 . 5 ) 









The background equations (3.4) for this Lagrangian are 


3Mp\H 2 = 2qP q — P + V , (3-6) 

Mp\(3H 2 - 2 HH') = -P + V , (3.7) 

- 3v^P, = • (3-8) 

The total derivative of the equation (3.7) with respect to IV becomes 

y' = + 6 qP q - 2M P \(H' 2 + HH") . (3.9) 


The quantities (2.8), (2.9) and the square of sound speed (2.10) 




77 _ m 2 _ 

J-s — M Pl e — -jr^, 

2 _ Pq 

° s ~ P q + 2 qP m ’ 


Qs = 


q(Pq + 2 qPqq) 

H 2 


(3.10) 


where we have used the equations (3.6) and (3.7) for Ps- Substituting them into the field 
equation (3.8), we obtain 


q'P q = P' = Mp, • -j-f-riH' 2 + HH") = l'{N) , 
c* + 1 

P(N) = I(N) + C , (3.11) 


where C is an integration constant. If we provide explicit functions for H(N) and c s (N), we 
can obtain the function I(N) as an explicit function of N after integrating (3.11). Substi¬ 
tuting (3.11) into the background equations (3.6) and (3.7), we get 


l 

Q 


I' 

Mp { HH' 


= AN) , 


In— = J(N) , 

qo 


(3.12) 


where qo is a constant which has the same mass scale as q. Then, if we can invert the function 
J(N) as N(q) = J^ 1 [ln(g/^ 0 )], we arrive to 


P(q) =I°J 1 [ln(g/? 0 )] + C , (3.13) 

V(N) = I(N) + C + Mp l (3H 2 (N) - 2HH’(N)) . (3.14) 


Here, we will set C = 0, since C always vanishes from the Lagrangian P — V. By plugging 
now q and N into X and ip and using the field redefinition (3.3), we finally obtain the form 
of the original Lagrangian, as 

p(x) - v^)) = ioj~ 1 [Hx/qo)} 

-I{p>m-M 2 l (3H 2 (v{ct>))-2HH^((t>))) . (3.15) 

The Lagrangian (3.15) created by using this formalism 2 will produce H(N) and c s (N) which 
we have provided at (3.11). If we tune H(N) and c s (N) so as to satisfy the observational 

2 A different formalism for the reconstruction of this Lagrangian can be found in [18, 19]. These authors 
use another affine parameter, instead of the e-folding number N. 


- 6 - 









constraints (2.18) and the stability condition (2.11), the reconstructed Lagrangian fulfills 
these constraints automatically. And if we use this formalism, we can ensure the stability of 
the system from the start, a very important property which cannot be generically ensured 
by using other methods. Hence, with our procedure it becomes relatively much easier to 
construct stable models than with other techniques. 

We shall now make a comment on the specific case V = 0. In this case, we obtain 


I' = - 2 H’ 2 - 2 HH”) 

from (3.14), and substituting this into (3.11), we get 

2 H’ H" 

° s ~ ~ 1 + 3H + 3H ' 


(3.16) 


(3.17) 


If V = 0, c s and H are not independent from each other and, therefore, we need just to 
provide either c s (N) or H(N ) in order to determine the unknown function P(X), using the 
same procedure as (3.12)-(3.13). This exactly corresponds to a previous result by [23]. 

3.1.1 Example 1 

We will now consider a specific model for the Lagrangian P(X ) — V(<f). Let us consider two 
explicit functions for P[(N ) and c s (lV), namely 


H{N ) = Ho(4N + p) 4 , 

T 


c s (N) = 


= const. 


(3.18) 

(3.19) 


where p and Hq are constants, to be determined by observations. The values of (e, f s , g s ) 
for this model are 

P 


e = 


fs — 9s — 


41V + p 

and the values of the observables (2.15)-(2.17), 

2p + 4 16 p 


4 N+p’ 


n, - 1 = - - 


41V + p 1 


r = 


VS(4 N+pY 


a.s = - 


8 p + 16 

(41V + pf 


(3.20) 


(3.21) 


The tensor to scalar ratio r is smaller by a factor y/3 than the value when cj. = 1. If we 
take N = 60 as the time when the curvature perturbation whose co-moving wave number 
is the pivot scale A;* = O.OSMpc' 1 grew up to the scale c s /.Hjjv = 6o, the observable (2.18) 
provides the following constraint on p: 1 < p < 2.5. Hereafter, we will always take N = 60 as 
the time when the curvature perturbation whose co-moving wave number is the pivot scale 
7’* = O.OSMpc” 1 grew up to the scale c s /H |a/=6o- The stability condition (2.11) becomes 


> 0 


F s = Ml —- 
A P1 4 N+p 

Qs~Ps= (4" 1 l^s = 2P S >0 . 


(3.22) 


For 0 < N < 60, p > 0 or p < —240 is required from (3.22). Combining these constraints, 
the parameter p must satisfy 

1 < p < 2.5, (3.23) 
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during inflation. On the other hand, Hq is also determined from the observational value of 
(2.14) coming from the Planck data, as 


Hq = Mpi 


2ny/2p 


__-p 

3V 4 (4-60 + p) E ¥ C 


1/2 


(3.24) 


Now, we reconstruct the Lagrangian using the two explicit functions (3.18) and (3.19). 
Using (3.11) and (3.12), we find 


P = ^-HH' 


(3.25) 


and 


In— = -In 


qo 


HH' 


ill! (p(4N + p) 


y -2 
2 


(3.26) 


We can fix the integration constant in J(N) to Hq without loss of generality. If we take 
p = 2, q becomes a constant and we cannot invert about N. For this model, nevertheless, we 
can reconstruct the kinetic function P(q) without inverting q in N(q) for all values of p, since 
both q and P are functions of HH'. Solving HH' about q and substituting it into (3.25), we 
obtain 

. . MtoFa / o \ 2 

p(q) = - • (3-27) 


i \qoj 

On the other hand, the potential term becomes 

V = \ml,Hl(AN + p) 7 i (l-|) 
The field redefinition (3.3) for this model yields 

2 p 1 / 4 ^ 


H 0 (p-6) 

and substituting this into (3.28), we obtain 


(4N + p)~— 


(3.28) 


(3.29) 


*+ ) = p {-W>VG* 


4 p 

p—6 


1 - 


e( 0 ) 


(3.30) 


Inserting q into X and fixing the integration constant as qo = Mp\HqM 2 /\/2, where M is a 
constant, we finally obtain the reconstructed Lagrangian 


p(x) - v{4>) 


X 2 


f p — 6 


—j - 3 AH, Hi -2--6 

M 4 P1 V\/2p 1 / 4 AfpiM 


4 p 

p—6 





(3.31) 


This model does not have a canonical kinetic term since we have provided the constant sound 
speed c s ^ 1, and the potential term is a power of the scalar field and a monomial, if we 
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ignore the slow-roll parameter in the potential. Moreover, from (3.26), the kinetic term is 
proportional to H 2 e. As in (3.28), the potential term is roughly proportional to the Hubble 
parameter; V ~ 3Mp^H 2 and, therefore, this model inflates the universe through slow-rolling 
of the scalar field, as expected from the form of the Lagrangian (3.5). We note that the 
constant M does not affect the evolution of the system since we can remove it from the 
Lagrangian by defining a dimensionless scalar field as (f> = (j>/M. We need some mechanism 
for the reheating process after inflation ends. There are some works about reheating with non- 
canonical kinetic terms by gravitational particle production (see, e.g., [26]), but a detailed 
analysis of the reheating stage is beyond the scope of the present paper. 

If, in particular, we take p = 2 and ignore the slow roll parameter e in the potential, we 
obtain a quadratic potential 

yquadra (^) _ 12AIp]H^ 

The value p = 2 is consistent with the astronomical observations and with the stability 
condition, and hence this inflation model with quadratic potential can survive if we modify 
the kinetic term of the model properly, such as we have done. Of course, the model we have 
considered here is just an example and we can easily construct other models in the same fash¬ 
ion which are stable and consistent with observations by using this powerful reconstruction 
formalism. 


(3.32) 


3.2 G-inflation 

In this subsection, we consider a different Lagrangian, namely 

K(<f>,X) = P(X) , G(X) = G{X) . (3.33) 

The background equations (3.4) for this Lagrangian become 

3Mp\H 2 = 2 qP q -P + m^2q*G q , (3.34) 

M P i(3H 2 - 2 HH') = -P- H^lqqG q , (3.35) 

(\/2 ~qPq + 6 HqGq)' — 3{\f2qP q + 6 HqG q ) = 0 (3.36) 

and this field equation tells us that 

y/YqP q + 6 HqG q = De 3N , (3.37) 


where D is an integration constant which has mass dimension 2. If q has a constant value as an 
initial condition, the Lagrangian (3.33) is also constant and the universe evolves as a de Sitter 
spacetime, H = const. In such situation, the constant D must be 0 because the derivative 
term in the field equation (3.36) vanishes. In what follows we will discuss reconstruction for 
the case when q does not have a constant value and D ^ 0. For a formulation of the case 
D = 0, see the Appendix. 

Using (3.37) and the background equations, we can write all quantities q(N ), H(N ), 
and their derivatives with respect to N as 


P = -3M^H 2 + D\/2qe 3N , 


q'Gq = a 


2 M 2 x H' De 3N 

VM^)~ H 


(3.38) 

(3.39) 
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and the quantities Ps and Gs are now 


Xs = P s (q, g', g"; N), Gs = Gs(q , g'; N) . (3.40) 

If we provide explicit functions for H(N ) and c s = yj 'Ps/Gs, as in the case of P(X) — 
we have here to solve a nonlinear second order differential equation on q(N), i.e. 
c^(N) = P s (q, q'; g"; N)/Gs(q, ( l\ X). Then, after obtaining two independent solutions of 
the differential equation, say gi(IV) and q2(N), we invert them as N\(q) and A^(g). And 
substituting either of them, N\(q) or A^g), into (3.45) and (3.46), after integration we will 
get the explicit functions P(q) and G(q). Inserting now q into X , we will finally obtain the 
so-called reconstructed Lagrangian. 

However, two difficulties can be envisaged in this procedure, in spite of the fact that 
one can perform it step by step, in principle. First, it may generally be difficult to solve 
the nonlinear second order differential equation analytically to find the two solutions explic¬ 
itly, and hence we may not be able to obtain analytic functions for the Lagrangian to be 
reconstructed. Second, in the case when we do have a explicit function for c s (N), we can 
ensure the subluminality of the propagation speed of the scalar field, but it is generically 
very unclear whether the stability condition (2.11) is satisfied or not. To confirm stability, 
we need to obtain rather complex forms for Pg and Gs using the solution qi(N) or g 2 (-/V) 
after solving the differential equation on q(N). Therefore, this is usually not a good way at 
all to provide a explicit function for c s (N) directly. 

To confirm stability, it is more reasonable to tune Ps and Gs as they satisfy the stability 
condition (2.11), instead of giving a explicit function for c s (N) directly. In order to impose 
the stability condition (2.11) on models in a more convenient way, we first derive a relation 
between H 1 and a product made of the undetermined function G(q), instead of c s (N), as 

M^H' = -y/2q*G q . (3.41) 


In this case, the expression (3.46) turns into a first order differential equation on q(N), namely 


D y/%^ N = M^HH' 



(3.42) 


If we now provide a explicit function for H(N ), we have more chances to solve this differential 
equation without problems. This equation is still a nonlinear equation, but now it is first 
order, after setting the relation (3.41), and thus much easier to solve, in general, than the 
second order one. Furthermore, the quantities Ps and Gs become 


Ps 


(1 + e) 2 ’ 


Gs 


3 Mpje 
1 + e 




(3.43) 


where we should use the solution of (3.42) for q/q' in Gs- These expressions for Ps and Gs 
are relatively concise and thus we can constrain models to satisfy the stability condition quite 
easily. The stability condition (2.11) yields in our case 


e' > 0, e(l + e) f 2^ + 1 j - e' > 0 . (3.44) 

After we find a proper function for H(N) which satisfies (3.44) and the observational con¬ 
straints (2.18), we invert q(N) as N(q), substitute it into (3.45) and (3.46) after integration 
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on N, and insert back q into X, to finally obtain the reconstructed Lagrangian: 


C = P(X) - G(X)g^X, l V^ , 

P(X) = -3M^H 2 (N(X)) + DV2Xe 3N( ' X) , 



2 M^H' 



N=N(X) 


(3.45) 

(3.46) 


It is guaranteed that from the Lagrangian created in this way, the function H(N) we had 
originally provided in (3.42) will be derived, and thus the resulting system will be stable. We 
note that the formalism we have discussed here is just one way to construct a stable system. 
We could assign a different relation related c s , instead of (3.41), to determine the unknown 
functions P(X) and G(X). The reason we use the relation (3.41) is because it results in the 
tractable expressions (3.44) for the stability condition. Below, we will construct a specific 
model of G-inflation using the above formalism. 


3.2.1 Example 2 

As a concrete model of reconstruction of G-inflation, we consider the following one 

H 2 (N) = HqQ 2iN + A 2 , (3.47) 


where Hq and A are constant with dimensions of mass, while 7 is a dimensionless constant. 
We can solve the equation (3.42) for this model and the solution is 


Here, 


q(N) = Z) 2 e (4T 


6 )JV 


D 2 = 


8A/ 4 ,iL 4 

D 2 


7 2 ( 7 -l ) 2 . 


e 


7 

\ + x' 


M 2 

= H^ N 


j = 2r r 2 x 
(l+x) 2 ’ 

>0 , 


and the stability condition reads 


(3.48) 


(3.49) 


2^2 (7 + ix + 7) (^ 4 ) - 27 2 x 

, 7 > 0 , -- > 0 

(l + 7 + x) 2 (l+x + 7) 2 


(3.50) 


If 


2 -y /2 3 2 + V 2 

0 < 7 < —g— ° r 2 “ 7 - —2— 


(3.51) 


the conditions above are always satisfied and the model is stable. On the other hand, if we 
restrict to values of the period id, the parameter region 7 < — 1 is also stable, although 
in this case the universe is in the phantom phase, where the null energy condition (NEC) is 
violated: H' < 0 (H > 0). However, x should grow exponentially when N increases provided 
7 is negative, and it would become much larger than 1 at very early times. Therefore, we 
have to set the initial condition such that x < 1 in order to obtain a stable NEC violating 
universe. 
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(3.52) 


The values of (/ s , g s ) for this model are 


fs = 


27 s + 27 — 27X 


9 s = 


27 X 


1 + 7 + x ’ 1 + 7 + x’ 

and the values of the observables (2.15)-(2.17), 

27 37 2 + 37 — 'yx 


n, - 1 = —: 


r = 


1 + x 1 + 7 + x 
327 2 x 


/ ( 27 s — 3y)x 
(1 + 7 + x) 5 / 2 V 3(7 — 1 ) 


OLg — 


27 127 2 x 


(1 + x) 2 (1 + 7 + x) 2 


(3.53) 

(3.54) 

(3.55) 


For simplicity, we shall consider the cases x <C 1 or x 1 at IV = 60 only. If we consider the 
case x 1 at IV = 60 and discard all the other terms without x, the spectral index becomes 


n s - 1 ~ 7 . 

From the observational data (2.18) we have that 7 ~ —0.032. However, this value of 7 is 
excluded by the stability condition (3.51), and thus we cannot obtain any stable inflationary 
model in the region i»l with the Hubble parameter (3.47) and the relation (3.41). On 
the other hand, if x <C 1 at N = 60 and we ignore the terms proportional to x, the spectral 
index (3.53) becomes 

n s - lr - 57 , 

and the observational value for spectral index in (2.18) constrains 7 to be 

0.00494 < 7 < 0.00746 . (3.56) 

This parameter region of 7 can indeed satisfy the stability condition (3.51) and also other 
observational constraints in (2.18), as 

r ~ 32 7 5/ V /2 < 0.0002 , 
a s ~ —2y, —0.01492 < a s < —0.00988 , 


Thus, we find values for the parameter 7 consistent with the observations and ensuring the 
stability of the model when id. The tensor to scalar ratio is very much suppressed and 
this model has little amount of primordial gravitational waves in the consistent region of 7 , 
(3.56). 

Hereafter we will consider the case i< 1 and 7 = 0.006 only. It is important to note 
that 7 is determined to be positive and, hence, the NEC violating universe (in which 7 has 
a negative value) is excluded by the Planck observations. 

In the case i< 1, the Hubble parameter is dominated by the first term on the right 
hand side of (3.47) and that term is responsible for the energy of inflation. Thus, we can 
evaluate the values of Ho and A using the observational value for the scalar spectrum from 
Planck, as 


H 2 e 2 ^ N 

TV ~ —-- 

C 87 t 2 c s J 7 5 



JV=60 

2.1 x 10 25 GeV 2 > H 2 , 


(3.57) 
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where we have used the approximations J~s 2M| 1 7 2 x and Qs ~ 2Mpj7- (3.57) yields 
Hq <C 10 13 GeV, and the energy scale of inflation is relatively smaller than the typical energy 
scale for the potential driven inflation, as in Example 1. If we set Hq = 10 10 GeV, we can 
determine A as A = 1.7 x 10 8 GeV from (3.57). 

Now, we reconstruct the Lagrangian by substituting (3.47) and (3.48) and its inversion 
N = N(q) into (3.45) and (3.46), with the result 

C = P(X ) - G(X)g^W^^ , 

P(X) = MpjTfg (4 7 2 - 47 - 3 ){X/D 2 )H^ - 3M^A 2 , 

G(X) = ~ 2 ^ M ^ Ho ^/(x/D 2 )^ + A2/ff2 

A 2 \”/24 9 \(W» 2 )^+A 2 /<T 
' P V Hi) \ n ) 499- 2n 

for 7 = 0.006. Note that the quantity D does never affect the evolution of the system since 
we can remove it from the Lagrangian by rescaling the scalar field as (f> = (j>/D , similarly 
as in Example 1. The Lagrangian obtained here does not include a canonical kinetic term 
and has a very complicated form. In general, the reconstructed Lagrangian does not always 
show the familiar form which includes a canonical kinetic term plus irrelevant terms. How¬ 
ever, this Lagrangian will certainly produce observables which are fully consistent with the 
observational values and, moreover, the cosmological model will be stable. 

We should stress again that ensuring the consistency with observations and an easy 
determination of the stability of the model are two of the main advantages in favor of using 
the reconstruction formalism over other methods. In exchange for these consistency and 
stability, one has to be ready to make sometimes sacrifices in terms of the simplicity of the 
Lagrangian. However, taking the lesson from effective field theories, the complexity of our 
Lagrangian can be explained in terms of the existence of an unknown (more fundamental) high 
energy theory which becomes valid at an energy scale much higher than the one corresponding 
to our ‘effective field theory’, and we do not have to worry so much about the complexity of 
our Lagrangian, necessarily. Again, the model we have shown here is just an example, and 
we are in principle able to reconstruct a Lagrangian which has a canonical kinetic term by 
using a different H(N). The important thing to be noted is that we have presented in detail 
a full-fledged formalism of reconstruction which can be applied to any function of H(N). 

4 Conclusions 

We have developed in this paper a reconstruction formalism for two different types of La¬ 
grangian, namely P(X) — V(<f>) and G-inflation, with which we can make the models con¬ 
sistent with the observational data from the Planck mission and evaluate the stability of 
the models easily. We have reconstructed two specific examples for each type of Lagrangian 
which are safe from any instability and satisfy the observational constraints for the spectral 
index, the tensor to scalar ratio, and for the running of the spectral index. For the case 
of P(X) — V(4>), although it exhibits a quadratic potential for q i, we were able to obtain a 
model consistent with the observations by adjusting the sound speed c s to be smaller than 1 , 
namely c s = yH/3. On the other hand, for the case of the G-inflation, we have been able to 


(3.58) 


(3.59) 
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construct a stable model in a way which is much easier than with other methods, all thanks 


to the reconstruction formalism we have developed in the paper. 

We insist that the reconstruction formalism we have put up here is a very general 
procedure which can be applied, in principle, to any model if both H(N ) and a quantity 
related to c s (N ) are provided. It is well known by now that using this and other existing 
formalisms, one always ends up with a Lagrangian which is consistent with the astronomical 
observations when one inputs the explicit functions of H(N) and of the quantity related 
to c s (N) from the observational values of the power spectrum, (2.13) and (2.14), and so 
on. However, after applying these methods the Lagrangian one obtains does not necessarily 
satisfy the stability condition, a condition that is generically quite difficult to impose. And 
if the Lagrangian does not satisfy the stability condition, this means that the theory cannot 
actually interpret the observations, and we can exclude it. In our case we impose this 
condition at an earlier stage of the procedure, so that it is guaranteed in the end results and, 
hence, we can judge easily whether the resulting theory is appropriate to describe inflation 
or not. Further, our formalism does not make use of the slow roll approximation, and for 
this reason it can also be employed with all families of models in which the scalar held does 
not evolve with slow rolling. 

Finally, in this paper we have not considered the reheating process and we have also 
not included ordinary matter in the action. To include matter is quite important, not only 
for dealing with the reheating process but also in order to be able to construct dark energy 
models using this reconstruction formalism. We leave these issues for a future work. 
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(Spain), projects FIS2010-15640 and FIS2013-44881, and by the CPAN Consolider Ingenio 
Project. 

A Reconstruction of the exact de Sitter solution in G-inflation 

In the case of exact de Sitter background and constant velocity of the scalar held, i.e. H = 
const., q = const., the parameter D must be zero for consistent evolution, and all background 
equations become algebraic equations of q. At this time, the background equations (3.4) 


become 


3M^H 2 = -K , 
y%K q + 6H^2^G q = 0 . 


(A.l) 

(A.2) 


We obtain a variation of H(q) by q from (A.2) as 


SH _ K qq Kq K q G qq 


(A.3) 


5q 3y/2qG q 6-^/2 q 3 G q 3y/2qG^ 


Using this variation, the coefficients of the linear scalar perturbation become 



(A.4) 
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This expression is the generalization of the corresponding expression in [15]. In order to obtain 
a reconstructed Lagrangian, we have to provide two inputs, as the functions of q, H(q ) and 
(? s {q), instead of H(N ), and a quantity related with c s (N), as the stability condition (2.11) 
is satisfied. For example, if we provide 


H(q) = 


cl = 


/ 1 1 

( q ^ 

/ 3Mpj 1 

l 2 Mj s ,d 

y(i - y) 


6(1 - y /2){l-y+ 1/(1 -y/2)} 

q 


Ml 


dsV 


0 < y < 1 


(A.5) 

(A.6) 


where M^s and y are constants which have a mass dimension. The function K(X) is trivially 
reconstructed from (A.l), as 


K(X) = -X + 


X 2 


2Mj S n 

and we obtain the quantities J~s and Qs using K(X = q), as 


(A.7) 




Gs 


Mp { g y(l-y) 

0 2 6(1 -y/2) ’ 


(i-y + V(i - y/2) - 2 q^j 


(A.8) 

(A.9) 


Comparing (A.8) and (A.9) with (A.6), we find G qq = 0, and hence the reconstructed function 
of G is 


G(X) = DX , 


(A.10) 


where D is an integration constant. We ignore the second integration constant since eventu¬ 
ally it becomes a total derivative term in the action. If we set D = 1/M% S , the background 
equation (A.2) determines the value of //, 


6// 2 (1 - y c ? 

M 2 ! y 2 (l-y c /2)’ 


(A-11) 


where y c is just a number which has the same value as y, but its variation with respect to q 
is zero, 5y c = 0, while 8y = 5q/M^ s fi. This result is completely consistent with [15]. 
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